A simple mean field equation for condensates in the BEC-BCS crossover regime 
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We present a mean field approach based on pairs of fermionic atoms to describe condensates in the 
BEC-BCS crossover regime. By introducing an effective potential, the mean field equation allows 
us to calculate the chemical potential, the equation of states and the atomic correlation function. 
The results agree surprisingly well with recent quantum Monte Carlo calculations. We show that 
the smooth crossover from the bosonic mean field repulsion between molecules to the Fermi pressure 
among atoms is associated with the evolution of the atomic correlation function. 
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Recent studies on ultracold Fermi gases and molecular 
condensates Q address an intriguing topic, the crossover 
from a Bosc-Einstein condensate of composite bosons to 
a fermionic Bardeen-Cooper-Schrieffer superfluid (BEC- 
BCS crossover) Q. By magnetically tuning the interac- 
tion strength near a Feshbach resonance [3] , a molecular 
BEC can be smoothly converted into a degenerate Fermi 
gas and vice versa. Experimental 0,01,01 an d theoretical 
research 0, II] into the quantum gases in the crossover 
regime are highly active and may provide new insights 
into other strongly interacting Fermi systems. 

In contrast to weakly interacting atomic BECs, for 
which a simple mean field description based on the Gross- 
Pitaevskii equation has been very successful , theoreti- 
cal models on the condensates in the crossover regime are 
generally very sophisticated and require expertise bor- 
rowed from condensed matter theory. The difficulty in 
providing a simple model for the fermionic system comes 
from, firstly, the lack of a small expansion parameter. 
The full range of atomic scattering length a should be 
taken into account to describe the crossover. Secondly, 
quantum many-body correlations are intrinsically more 
complicated for fermionic systems than for bosonic ones. 

The BEC-BCS crossover, however, suggests an alterna- 
tive approach to model the strongly interacting fcrmions 
based on composite bosons. This is possible since a Fermi 
gas in the crossover regime constitutes the same quantum 
phase as of a condensate of interacting pairs. Recent ex- 
periments on the wave function projection H and on the 
pairing gap 01 indeed indicate that near the Feshbach 
resonance, a large fraction of fermionic atoms are paired 
at low temperatures. From these observations, we pro- 
pose a bosonic mean field equation, complementary to 
the fermion-based BCS approaches, to describe the atom 
pairs in the crossover regime. Our mean-field approach 
is relatively simple and well-behaved near the resonance. 
We obtain analytic expressions for the chemical poten- 
tial and the equation of states, which agree very well 
with other calculations. In particular, we find the chemi- 
cal potential in the unitarity limit is ~ 0.4357 times that 
in the BCS limit, in excellent agreement with the recent 
quantum Monte Carlo calculations of 0.42 ~ 0.44 |T3.[T^. 

We consider an ultracold gas of two-component 
fermionic atoms. At low temperatures, only atoms in 



different internal states can pair via s-wave interaction. 
For simplicity, we assume the interaction range is zero. 
In the absence of many-body effects, the center-of-mass 
motion of an atom pair ^Iq(R) is decoupled from the in- 
ternal relative atomic motion i/jo(r) = (47rr 2 ) _1 / 2 -0o(r) 
with r = \f] the atomic separation. Given the atomic 
scattering length a, ipo (r) satisfies Schrodinger's equa- 
tion, 



-<(r) = -E b Mr) 



(1) 



with the boundary condition ipo(0) = —aip' (Q). Here m 
is the atomic mass, 2nh is Planck's constant, and Eb is 
the molecular binding energy 

For positive scattering lengths a > 0, the bound state 
is described by ip = (2/a) 1 / 2 e" r / a with E b = h 2 /(ma 2 ). 
The size of the molecule is given by (r) = a/2. For 
negative scattering lengths a < 0, the bound state does 
not exist and the ground state energy is — Eb = 0. 

Now consider a condensate of pairs with a density dis- 
tribution n(R) in a slow-varying potential well V(R). We 
introduce the many-body wave function to include the 
condensate of the bosonic pairs $(R) = n(i?) 1 / 2 as well 
as the internal atomic correlation - 0( r )- The mean field 
equation for the composite bosons is then 



4m 



h 2 df. 



\-V + U)V(R)ip(r) 
ip{0) = -ad r i>(0). 



,*(Ji)V(r)(2) 
(3) 



Here \x m is the chemical potential, U = g\^(R)\ 2 is the 
mean field interaction and g is the interaction term. 

In conventional approaches, g is given by the scatter- 
ing length of the bosons. For pairs of fermions, scatter- 
ing length a m is determined by that of the constituent 
atoms as a m = 0.60a, resulting from an effective repul- 
sive potential between molecules ^4|. This dependence 
can be understood in a simple picture. Low-energy col- 
lision with a repulsive interactions acquires a scattering 
length which is proportional to the size of the scatterer. 
For pairs of atoms, we have a rn ~ (r) = a/2. 

From the above considerations, we hypothesize that 
the interaction term g is effectively proportional to the 



2 



interatomic separation r as 

H 2 

9 = g( r ) = c—r, (4) 
m 

where c is a dimcnsionless constant. 

To proceed with minimum algebra, we consider a uni- 
form gas with a density \^(R)\ 2 = n. Eq. (2) becomes 

(-^d 2 r +gW(r)=ii m i>(r). (5) 
m 

To determine c, we consider the BEC limit (na^ <C 1), 
where the mean field term can be treated perturbativcly. 
That is, the expectation value of U based on the bare 
molecular wave function ip (r) = (2/a) 1 / 2 e" r / a should 
yield the molecular mean field shift 4irh 2 a m n/2m, 



since the pair is compressed. This result also provides 
a simple picture to understand the augmentation of the 
molecular binding energy reported in Ref. 0| ■ 

We extend the mean field model to the crossover and 
the BCS regime, where the atom pairs strongly overlap. 
Although it becomes less clear if the mean field approach 
can fully capture the Fermionic nature of the gas, our aim 
here is to determine an effective potential which can best 
describe the system in the BEC-BCS crossover regime. 

In this regime, the four-body calculation of a m = 0.6a 
is no longer valid, and we determine the mean field in- 
teraction U from the properties of the Fermi gas. First 
of all, in the dilute gas limit, we still expect the inter- 
action to be proportional to the square of the bosonic 
field, U — g\^(R)\ 2 . Secondly, we exploit the asymptotic 
behavior of the gas in the weak coupling limit na 3 — > _ , 
where the system approaches an ideal degenerate Fermi 
gas with the chemical potential 



2i:h 2 a m n 



(6) 



Using a m = 0.60a, we find Eq. (6) can indeed hold for 
arbitrary scattering lengths a m based on the linear mean 
field potential in Eq. (4). We determine c = 47ra rn /a « 
7.5. 

From Eq.(3), (4), and (5), the exact solution of the 
pair wave function is given by 



NAi{c 1/3 n 1/3 r 
—a d r ip(0), 



C-' 2/3 LL 



,/Eo) 



(7) 

(8) 



where N is the normalization constant, Ai(x) is Airy's Ai 
function, and Eq = h 2 n 2 ^ 3 /m. Notice that the chemical 
potential (j, m in Eq. (7) is determined from Eq. (8). 



In the weak interaction limit < 



<C 1, the wave 



function ip(r) obtained from Eq.(7) is identical to the 
unperturbed one tpo(r) f° r r "C n^ 1 /' 3 . For r n^ 1 / 3 , 
ip(r) is exponentially smaller than ^(r) and approaches 



r^ 1 / 4 exp(— |r 3 / 2 ). This suppression for large atomic 
separation is expected since the interaction energy in- 
creases when the pairs start overlapping. As a conse- 
quence, the pair wave function ip{ r ) is compressed to a 
smaller size than that of a bare molecule. Similar effect 
is also discussed in Ref. [THj 

The distortion of the pair wave function can be char- 
acterized by an effective shift in the binding energy Ef,. 
In the weak interaction limit, the shift can be defined as 



lim fx m = 2E F 



(6n 2 n) 2 / 3 h 2 



(10) 



where E F = h 2 k F /2m is the Fermi energy and k-p = 
(S^n) 1 / 3 is the Fer mi wave number. 

Based on Eq. (5), we find that the above density de- 
pendence fj, m oc r?l 3 can be satisfied only when the in- 
teraction term g is again a linear function of r. Taking 
the limit of a = 0~ and assuming g(r) = c'(h 2 /m)r, 
we can solve the chemical potential from Eq.(7) and 
Eq.(8) as [i m = ac' 2 / 3 E , where -a ps -2.338 is the 
first zero of the Ai(:r) function. Equating \i m to 2E-p 
yields c' = 6ir 2 a~ 3 ^ 2 « 16.56. This value is about twice 
as large as c. 

We first test the equation in the unitarity limit a = 
±oo. Fermi gases in this limit have been extensively stud- 
ied, for which a universal and fermionic behavior is ex- 
pected . Due to the divergence of the scattering length, 
we expect the only energy scale in the system is the Fermi 
energy Ef. From the boundary condition d r ip(0) = and 
g = c' (h 2 /m)r, we determine the chemical potential as 
fx m = a'c' 2 / 3 Eo, where —a' « —1.019 is the first zero 
of the Ai'(x) function. Given c' = 67r 2 a -3 / 2 , we get 
pL m /2 = (a /a)E F w 0.4357-Ef- This result agrees ex- 
cellently with recent quantum Monte Carlo calculations 
which gives fx m /2 = 0.U(1)E F ^ and 0.42(1)£ F []]|, 
and the measurements 0, where the uncertainties 
are larger. We, however, cannot exclude this agreement 
is coincidental. Near the unitarity limit, we have 



r (r)(- 



h 2 dl 



)i){r)dr = -E b (l + o{na 3 )). (9) 



The binding energy correction o(na 3 ) is positive. This 
effect is absent in the calculations for point-like bosons 
|l6j since it originates comes from the internal degree 
of freedom. This increase in binding energy is expected 



2E F 



a 
a 



a 



-1/2 



a'k F a 



(11) 



Next, we investigate the BEC-BCS crossover regime. 
Rewritting Eq. (7) and (8) using c! = &iT 2 a~ 3 / 2 , we get 
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1/k F a 

FIG. 1: Chemical potential fi m in the crossover regime (solid 
line). For large l/k^a, the chemical potential jix m approaches 
the energy of the molecular state — Et (dotted line). The solid 
dots and the open dots show the Monte Carlo calculations 
from Ref. 01 an d Ref- 0], respectively. 
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FIG. 2: Exponent 7 for the equation of states from the mean 
field calculation (solid line) , the BCS calculation (dashed line) 
[Tgl and the fit to the quantum Monte Carlo calculation (open 
circles) [HIH. The unitarity and BCS limit 7 = 2/3 and 
the BEC limit 7=1 are shown in dotted lines. 



ip(r) = NAi{a- 1/2 k F r~a-^-) (12) 

2E F 

k F a Ai(-an m /2E F ) 

«V2 = Ai'(-an m /2E F y [ ' 

The chemical potential /j, m calculated from Eq.(13) is 
shown in Fig. ^ We see that \i m approaches 2E F in the 
BCS limit and —E^ in the BEC limit, as expected. In the 
crossover regime, the values agree well with the Monte 
Carlo calculation from 01 ■ We can also evaluate the 
equation of states fi m + E(, oc n 7 , where the exponent 7 
pl ays a crucial role in the collective excitation frequencies 
From Eq.(12) and Eq.(13), we obtain 

dln(n m + E b ) 

7 = j; (14) 

a inn 

= 2 (i 1 ^V7i 1 ^ 3/2 WMA fl5 ) 

" 3 V 1 + fx m J V 1 + fc|a 2 + 2E F / n m ) 

The exponent 7 (see Fig. shows the expected behav- 
ior: 7 = 1 in the BEC limit and 7 = 2/3 in the BCS and 
unitarity limits. In the range of 1 < k F a < 00, 7 shows 
a dramatic variation. In the following, we show that this 
dramatic variation is directly linked to the crossover na- 
ture of the quantum gas and is a result of the distortion 
of the pair wave function ip( r )- 

From Eq. (12), we calculate ijj(r) for k F a = 1/2, 2, and 
±00, shown in Fig. [3] For k F a = 1/2, we see very small 
deviation of ip(r) from the bare molecular wave function 
ipo(r). For k F a = 2, ip( r ) is clearly different from ipo(r) 
with a higher probability amplitude for r < k^ 1 and a 
lower amplitude for r > fcp 1 . This is the compression 
effect we discussed. In the unitarity limit k F a — ±00, 
the atomic pairing is fermionic since bare molecules dis- 
sociate at this point. The mean atomic separation of 




r (k F 1 ) 

FIG. 3: Pair wave functions in the crossover regime. Wave 
functions il>(r) at fea = 1/2, kFa = 2 and kpa — ±00, shown 
in solid lines, are calculated based on Eq.(12). In the former 
two cases, the bare molecular wave functions ipo(r) (dotted 
lines) are shown for comparison. 

(r) w 2/fcp ~ O.bn" 1 / 3 suggests the size of the pairs is 
about half of the mean molecular spacing. 

The distortion of the wave functions leads to significant 
consequences for the quantum gas. Given the mean field 
energy as (U) oc n(r), the evolution of the pair size from 
(r) oc a in the BEC regime to (r) oc n -1 / 3 in the unitarity 
limit underlies the crossover nature of the interactions 
from the bosonic mean field repulsion (U) oc na to the 
Fermi pressure (U) oc ?i 2 / 3 . This explains the variation 
of the exponent 7 in Fig. (2). From these observations, 
we can qualitatively define the BCS regime to be k F a < 
0, crossover regime 1/2 < k F a < 00, and BEC regime 
< k F a < 1/2. The use of c' w 16.56 in the mean field 
term is appropriate in the BCS and crossover regimes 
and c w 7.5 in the BEC regime. 

The pair wave functions can be directly probed exper- 
imentally by radio-frequency (rf) excitations as demon- 
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FIG. 4: Bound-free Franck-Condon factors Ff(K) of the pairs 
for (from bottom to top) fea = 1/2, feo = 1, fcFffl = 2, 
ipa = ±oo (unitarity limit) and kpa — CP (BCS limit, dotted 
line). The arrows mark the peak positions K p k- In the inset, 
K p k is plotted as a function of 1/kva (solid line) together with 
the K p k for bare molecules (dashed line). 



stratcd in Rcf. |10ll22| . In these experiments, rf photons 
excite the bound pairs into another spin state in which 
no bound state exists. The excited pairs then dissociate 
into free atoms. Theoretical calculation based on bare 
molecules show that the excitation rate constant, or the 
bound-free Franck-Condon factor Ff (K) , reflects the pair 
wave function in the momentum space [2l|, 



Ff(K) 



irh 2 k 



sin (kr + 5)ipQ{r)dr 



, (16) 



where K = nrk jm, k, and d are the energy, relative wave 



number and the scattering phase shift of the outgoing 
atoms, respectively. 

To calculate Franck-Condon factors in the crossover 
regime, we replace ipo( r ) by ?K r ) an d assume the atoms in 
the outgoing channel do not interact 5 = 0. In Fig.0] we 
show that the Franck-Condon factors display a resonance 
structure in the crossover regime. The location of the 
peak Franck-Condon factor K p k provides a sensitive mea- 
sure of the atomic correlation leng th. In the BEC regime, 
Kpk approaches ^E b ^> E F |2lJ and suggests that the 
atomic separation is small compared to the intermolec- 
ular distance. In the crossover regime, K p k approaches 
a small fraction of Ep. The persistence of the resonance 
structure at unitarity and in the BCS regime indicates 
the correlation of the atoms in momentum space. This 
dependence is recently reported in [ToL ll ll| . A quantita- 
tive comparison with the measurements, however, must 
include the effects of the trapping potential and the fi- 
nite temperature [2^, which is outside the scope of this 
paper. 

In conclusion, we provide a simple mean field model 
to describe the BEC-BCS crossover. From this equation, 
many properties of the strongly interacting gas can be an- 
alytically calculated with high accuracy. Our model can 
also be easily generalized to include the external poten- 
tial and to study crossover effects in systems with lower 
or higher dimension. 
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